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ABSTRACT 



Context. Most of the Planetary Nebulae (PNe) observed are not spherical. The loss of spherical symmetry occurs somewhere between 
the Asymptotic Giant Branch (AGB) phase and the PNe phase. The cause of this change of morphology is not yet well understood, 
but magnetic fields are one of the possible agents. The origin of the magnetic field remains to be determined, and potentially requires 
the presence of a massive companion to the AGB star. Therefore, further detections of the magnetic field around evolved stars, and in 
particular those thought to be part of a binary system, are crucial to improve our understanding of the origin and role of magnetism 
during the late stages of stellar evolution. One such binary is the pre-PN OH231.8+4.2, around which a magnetic field has previously 
been detected in the OH maser region of the outer circumstellar envelope. 

Aims. We aim to detect and infer the properties of the magnetic field of the pre-PN OH231.8+4.2 in the H2O maser region that probes 
the region close to the central star. This source is a confirmed binary with collimated outflows and an envelope containing several 
maser species. 

Methods. In this work we observed the 61,6-52,3 H 2 maser rotational transition to determine its linear and circular polarization. As 
a result of Zeeman splitting, the properties of the magnetic field can be derived from maser polarization analysis. The H 2 maser 
emissions of OH23 1.8+4.2 are located within the inner regions of the source (at a few tens of AU). 

Results. We detected 30 H 2 maser features around OH231.8+4.2. The masers occur in two distinct regions that are moving apart 
with a velocity on the sky of 2.3 mas/year. Taking into account the inclination angle of the source with the line of sight, this corresponds 
to an average separation velocity of 21 km/s. Based on the velocity gradient of the maser emission, the masers appear to be dragged 
along the direction of the nebula jet. Linear polarization is present in three of the features, and circular polarization is detected in the 
two brightest features. The circular polarization results imply a magnetic field strength of |B|||~45 mG. 

Conclusions. We confirm the presence of a magnetic field around OH23 1.8+4.2, and report the first measurements of its strength 
within a few tens of AU of the stellar pair. Assuming a toroidal magnetic field, this imples B~2.5 G on the stellar surface. The 
morphology of the field is not yet determined, but the high scatter found in the directions of the linear polarization vectors could 
indicate that the masers occur near the tangent points of a toroidal field. 

Key words, masers, polarization, magnetic field, Stars: AGB and post-AGB 



1. Introduction 

During the final stages of their evolution, low and intermedi- 
ate mass stars (0.8-8 M ) evolve from the Asymptotic Giant 
Branch (AGB) phase to Planetary Nebulae (PNe). During this 
transition most of these objects lose their spherical symmetry. 
The process responsi ble for the change of m orphology is, so far, 
not well understood. Buiarrabal et al. (2001) have shown that, 
in most cases, the radiation pressure does not have enough en- 
ergy to drive the acceleration of the fast bipolar flows observed 
in pre-PNe (pPNe). Moreover, the mass ejection on a preferred 
axis requires an extra agent to collimate the flow. The possible 
mechanisms that could both provide the missing energy to drive 
the fast flows and that could collimate it on a preferred axis are: 
(i) a companion to the star (either a binary stellar companion 
or massive planet) and its tidal forces, (ii) disk interaction and 



(iii) magnetic fields - or a combination of these dBalick & F rank 
2002| |Franke7al][2007; Nordha us et all 120071 and references 
therein). 

From magnet o -hydro dynamic (MHD) simulations, 
iGarcfa-Segura et all d 19991) concluded that magnetic fields 
can indeed have a pronounced effect on shaping stars beyond the 
AGB stage and that, working together with rotation, they can 
produce collimated bipolar nebulae and jets. Since then, several 
other MHD simulations have been designed to investigate how 
magnetic fie lds can shape a pPN and, consequently, a Planetary 
Nebula (e.g., Garcfa-Segura et alJ2005HGarcfa-Dfaz et al.f 2008: 
iDennis et al.H2009h . 

Observations of the magnetic fields around evolved stars, 
however, are still rare. Most current magnetic field measure- 
ments are performed u sing Zeeman spli t ting (IZeemanl [T897t) 
of maser lines (e.g., Vlemmin gs et alJ 1200 lL 120061) . Each 
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maser species requires different physical conditions to occur. 
Therefore, different masers originate in different locations of 
the studied object. The SiO masers are found close to the cen- 
tral star (CS), in regions with a temperature of ~1300 K. The 
OH masers are generally expected to occur further out, at a few 
hundred AU from the CS. For AGB stars, the H2O maser emis- 
sion is expected to be intermediate between these two regions; 
but for a pPN/PN it can be f ound at a similar or greater radius 
than the OH maser emission jCohenll987HReid & Mentenll997t 
IVlemmings & Diamondll2006l) . 

In this work, we aim to investigate the polarization of the 
H 2 maser emission toward the pPN OH23 1.8+4.2 (Rotten 
Egg Nebula; Calabash Nebula), and infer the prope rties of its 
magn etic field. Located at a distance of ~1540 pc dChoi et al.l 
1201 ll updated from private communication), the bipolar 
nebula OH231.8+4.2 contains a binary system in its core, 
where a Main Sequence type A star accompanies an evolved 
star - the Mira variable QX Pup (Sa nchez Contreras et al.l 
2004). The SiO, H 2 and OH maser emission toward the 
Rotten Egg Nebul a have been observed several times be- 
fore (e.g.. iMorris et al.l 119821; iGom ez & Rodriguez 2001; 
Sanchez Contreras et al.ll2002b iDesmurs et al.ll2007t lEtoka et al l 
2008, 2009). These prior works have shown that the spatial dis- 



tribution of both the SiO and OH maser emission lie perpen- 
dicular to the nebular symmetry axis, consistent with the pres- 
ence of a torus around the equator of the central star. However, 
the H2O masers (located at only ~40 AU from the CS) do 
not follow the same pattern, and are instead located toward the 
bipolar structur e. Few of the prior studies inclu ded polarization 
measurements. ISanchez Contreras et al. I (12001 note that their 
SiO p olarization maps are not reliable and lGomez & Rodri guez 
(1200 lb found no circular po larization in the OH maser emis- 
sion. lEtoka et al.ld2008l 12009). however, found both circular and 
linear polarization in the OH masers, and conclude that a well- 
organized magnetic field seems to be flaring out in the same di- 
rection as the outflow. 

This paper is structured as follows: Section 2 describes the 
observations, data reduction, and calibration, and Section 3 in- 
cludes a presentation of the results. The results are discussed in 
Section 4 and the concluding analysis is presented in Section 5. 



2. Observations and Data Reduction 

The observations of the pPN OH231.8+4.2 were carried out 
on March 1", 2009, using the NRACQ Very Long Baseline 
Array (VLBA), under project code BV067A. We observed the 
H2O 6i > 6-52,3 rotational transition, at an adopted rest frequency 
of 22.235081 GHz. We used 2 baseband filters of 1 MHz 
bandwidth and, fo llowing the VLBA H2O observations of 
Desmurs et al. I (120071) (hereafter; D07), the filters were centered 
at Vi sr 44.0 km s _I and 26.0 km s _1 respectively. This was done 
in order to detect both the emission from the northern (NReg; 
~30 mas from the CS) and southern regions (SReg; ~40 mas 
from the CS) of the source. The observations were centered on 
RA (J2000) = 07 /, 42 m 16.93 and Dec (J2000) = -14°42'50".2. 
Both a lower and a higher spectral resolution correlation were 
performed. The former was undertaken in full polarization mode 
(to obtain all 4 Stokes correlations: RR, LL, RL and LR) over 128 
channels, with a resulting nominal channel width of 0.104 km/s. 
The latter correlation was performed in dual polarization mode 



1 The National Radio Astronomy Observatory (NRAO) is a facility of 
the National Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 



(to obtain RR and LL) over 512 channels, with a resulting nomi- 
nal channel width of 0.026 km/s. The synthesized beam size was 
~1.7x0.9 mas. The observations spanned 8 hours, of which 4.7 
hours were spent on OH231.8+4.2. 

For calibration purposes, we observed J0854+2006. This ref- 
erence source was used to perform bandpass, phase and polar- 
ization calibration. We used the Astronomical Image Processing 
Software Package (A IPS) and followed the data reduction proce- 
dure documented by Kemba ll et alj d 1995b . an d which has been 
successfully adopted by several authors (e.g., fviemming s et al.l 
l2001blSurcis etaLll201 lh . 

Using the low resolution data, we created image cubes for 
the Stokes parameters /, Q, U and V. The noise level measured 
from these image cubes ranges between ~3 and ~5 mJy/beam in 
the emission free channels. From the Stokes Q and U images, 
we derived image cubes for the linear polarization intensity (P) 
and the linear polarization direction (also referred to as electric 
vector position angle, EVPA). 

Using the high resolution data, we created image cubes for 
the Stokes parameters / and V. The noise level measured from 
these image cubes lay between ~6 and ~10 mJy/beam in the 
emission free channels. To estimate the polarization leakage 
and consequently the minimum believable fractional linear po- 
larization we determined the Stokes Q and U limits for the 
brightest maser feature with no detected polarization. We con- 
clude that, at 3cr, the minimum detected fractional polarization is 
~0. 1 %, a value that is consistent with the self-calibration method 
used during the data reduction, described in further detail by 
iLeppanen etail (1 19951) . 

In our analysis we adopted a signal-to-noise ratio cutoff of 
four times the measured noise (4cr). We define a maser feature 
to be successfully detected when maser spots located at similar 
spatial positions (within the beam size) appear and survive the 
threshold cut in at least 3 consecutive channels. The position of 
the maser feature is then define d by the emission in th e channel 
where the peak appears (see e.g. Richards et al. 20~PTt). 



3. Results 

3.1. Masers Location 

We detected 30 maser features in total, 20 of them located in 
the NReg and the remaining 10 located in the SReg. In FigureQ] 
we present four plots: l.a and l.b correspond to the NReg, and 
l.c and l.d to the SReg. Plots (a) and (c) show the distribution 
of all 252 maser spots that comprise the 30 maser features. The 
dashed vectors indicate the direction of the pPN jet. Plots (b) and 
(d) show the maser features, scaled by their flux densities (pro- 
portional to the size of the triangles), and color-coded according 
to their line-of-sight velocities. The black line segments, scaled 
in size by the degree of polarization, show the EVPA (see section 
13.21 ). The origin (0,0) of all four plots is centered on the brightest 
feature we detected, located in the SReg. 

In order to aid the comparison of the observations presented 
by D07 with the current work, we over-plotted both observa- 
tions in Figure [2] Since we did not use phase-referencing during 
our observations, unlike D07, we do not have a precise value 
for the absolute maser positions (but we do have accurate rela- 
tive astrometric information). Therefore, to be able to make this 
comparison, we assumed a common origin point at the center of 
the maser emission. We determined this center point by taking 
the mean position over all the 30 maser features detected here, 
and over 15 of the 16 features detected by D07. We chose not 
to include maser feature number 16 from D07 in the analysis, 



2 



Leal-Ferreira et al.: The magnetic field of the pPN OH231.8 




Right Ascension offset (mas) Righl Ascension offset (mas) 



Fig. 1. Plots a and b (top) corresponds to NReg, and plots c and d (bottom) to SReg. Plots a and c (left) show the position of all the 
measured maser spots, and plots b and d (right) present the maser features. On plots a and c the dashed vectors indicate the direction 
of the pPN jet. On plots b and d the size of the triangles are scaled by the maser fluxes, and the color scale is related to the velocity. 
The black lines represent the EVPA and intensity of the linear polarization. 



because it is clearly a spatial outlier in the field that contains the 
features we observed. Figure [2] shows the entire field of view 
(NReg and SReg). The filled triangles represent the 30 features 
detected in the current work and the hollow squares show the ob- 
servations from D07. Five regions (2.A, 2.B, 2.C, 2.D and 2.E) 
are drawn on the figure to facilitate their identification in the text. 

In the NReg, the positions of the brightest feature and its 
surrounding features, are similar to those found by D07, except 
that the current observations detected more widespread features 
(region 2.A). On the eastern side of the brightest northern fea- 
ture, at a separation of ~30 mas, we detected only one weak fea- 
ture, while two features were detected by D07 in this region, at 
~20 mas from the brightest maser (region 2.B). On the western 
side of the main feature we detected a group of 7 features that to- 
gether seem to compose a single extended structure (region 2.C). 
This group of features does not appear in the observations of 
D07. Instead, they found only two features in that same region. 



Moreover, the emissions we detected in the NReg are contained 
within a much wider North-South area (~26 mas) than that in 
the D07 image (~13 mas). 

In the SReg, the distribution of the features we observed is 
similar to those presented by D07, but appears to have moved 
away from the origin. Furthermore, we detected 10 features in 
this region, while D07 detected only six features (considering 
only regions 2.D and 2.E). The detections in the current obser- 
vations are spread over a wider East-West area (~36 mas) than 
that reported by D07 (~17 mas). 

The 30 maser features and their properties are listed in 
Table □ 

3.2. Linear Polarization 

We detected linear polarization for three of the 1 1 maser fea- 
tures in the SReg (S.05, S06 and S07). The degree of polar- 
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Table 1. Maser features detected on OH23 1.8+4.2 



Reg 


Label 


Aa (mas) 


AS (mas) 


/ (Jy/beam) 


.1 1 ( J y) 


V,,„- (km/s) 


P (%) 


EVPA (") 


NReg 


N.01 


-11.2057 


81.4210 


0.49 


0.81 


30.2 


— 


— 


NReg 


N.02 


+ 15.8905 


67.5530 


0.08 


0.12 


29.4 


— 


— 


NReg 


N.03 


-13.3862 


78.3560 


0.24 


0.70 


29.2 


— 


— 


NReg 


N.04 


-14.2320 


76.0570 


0.27 


0.32 


28.5 


— 


— 


NReg 


N.05 


-15.2988 


73.8400 


0.54 


1.08 


27.4 


— 


— 


NReg 


N.06 


-6.31373 


84.3500 


1.67 


1.61 


27.1 


— 


— 


NReg 


N.07 


-24.1890 


69.6120 


0.30 


0.46 


26.8 


— 


— 


NReg 


N.08 


-22.6036 


67.9370 


0.33 


1.33 


26.7 


— 


— 


NReg 


N.09 


-16.7119 


59.5600 


0.26 


0.38 


26.4 


— 


— 


NReg 


N.10 


-25.9197 


69.1120 


0.92 


1.72 


26.4 


— 


— 


NReg 


N.ll 


-21.0261 


65.8640 


1.42 


2.33 


26.4 


— 


— 


NReg 


N.12 


-8.18628 


80.2300 


0.72 


0.98 


26.4 


— 


— 


NReg 


N.13 


-15.4815 


58.4830 


0.34 


0.60 


26.2 


— 




NReg 


N.14 


-9.49584 


76.7490 


0.86 


1.34 


26.1 


— 


— 


NReg 


N.15 


-18.7994 


64.7640 


0.49 


0.79 


26.1 


— 


— 


NReg 


N.16 


-13.9092 


70.7840 


16.47 


19.24 


25.6 


— 


— 


NReg 


N.17 


-13.6490 


74.2530 


0.43 


0.73 


25.4 


— 


— 


NReg 


N.18 


-6.86105 


62.7080 


0.09 


0.11 


24.9 


— 


— 


NReg 


N.19 


-15.6860 


75.7780 


0.11 


0.13 


24.8 


— 


— 


NReg 


N.20 


-8.41686 


68.5520 


0.31 


0.40 


24.6 


— 


— 


SReg 


S.01 


-9.27134 


2.06101 


0.08 


0.11 


44.8 


— 


— 


SReg 


S.02 


-7.57631 


6.09103 


0.04 


0.04 


43.6 


— 


— 


SReg 


S.03 


-5.31655 


7.60502 


0.14 


0.23 


43.3 






SReg 


S.04 


-1.88907 


1.84101 


0.16 


0.31 


43.2 






SReg 


S.05 


-1.82468 


-0.05600 


2.55 


3.59 


42.5 


0.63 + 0.13 


-65 ±6 


SReg 


S.06 


-1.61411 


-2.82499 


1.35 


1.46 


41.9 


1.15 ±0.24 


+58 ±6 


SReg 


S.07 


+0.00000 


0.00000 


5.53 


6.77 


41.6 


0.28 ± 0.04 


-32 ±4 


SReg 


S.08 


-34.4941 


-1.73098 


0.10 


0.14 


40.2 






SReg 


S.09 


-29.9659 


-3.77899 


0.04 


0.05 


40.0 






SReg 


S.10 


-35.9028 


-3.57598 


0.05 


0.07 


39.3 







^,From left to right, the following maser parameters are presented: region of location (Reg); feature label 
(Label); right ascension offset, in reference to Figure[TJ(AQ'); declination offset, in reference to Figure[TJ(A<5); 
peak intensity (/); integrated intensity ( J I); radial velocity (V;, r ); P; EVPA. 



22.33 GHz 0H238.1 H 2 Filled Triangles: This work 

Hollow Squares: Desmurs et al. (2007) 




40 20 -20 -40 

Right Ascension offset (mas) 



ization and polarization angles of each feature are listed, with 
their respective errors, in columns 8 and 9 of Table Q] Those 
values correspond to the results given by the brightest channel 
of each feature. The P error is given by the rms taken from 
the P image, scaled by the intensity peak. The EVPA error was 
determined using the ex pression (Tevpa = 0.5 crp/P x 180°/7r 
dWardle & Kronbergl 19741) . 

In Figure[T]d, the black lines show the direction of the linear 
polarization of features S.05, S.06 and S.07. The length of the 
lines are proportional to the linear polarized intensity of each 
feature. 

In Figure [3] for each feature for which linear polarization 
was detected, we show the polarization vectors across all indi- 
vidual channels. Again the length of the vectors are scaled in 
proportion to the linearly-polarized intensity. 

3.3. Circular Polarization 

The total intensity (/) and circular polarization ( V) spectra of the 
maser feat ures were used to perfo rm the Zeeman analysis de- 
scribed by Vlemmings et al. (2002). In this approach, the frac- 
tion of circular polarization, Py, is given by 



Fig. 2. Overplot of the maser features observed by D07 and the 
features presented in the present work. The filled triangles show 
the 30 features we observed and the hollow squares represent the 
observations from D07. The dashed vector indicate the direction 
of the pPN jet. 



Pv — (Vmax Vmin) I I max 

= 2 xAf-f, x B^ Causs] /Av L[k „ 



(1) 



where V max and V m i„ are the maximum and minimum of the 
model fitted to the V spectrum and I max is the peak flux of the 
emission. Aj?_p is the Zeeman splitting coefficient, whose ex- 
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t'ea; arc S.05 



-1.4 -1.6 -1.8 -2 

Right Ascension offset (mas) 



22.23 GHz 011238.' 



1'cri.urc 3.06 




•1.4 -1.6 
Right Ascens 



22.23 GHz 0H238.1 1-1.0 



-1.8 

fset (mas) 




Feature S.07 




Fig. 4. Spectra of the Stokes / (top) and V (bottom). On the 
left hand side we present the spectra on the spatial position 
of the peak emission of the NReg. On the right hand side we 
show the spectra on the spatial position of the peak emission 
of the SReg. From the best fit model (blue line) we obtained 
B,|(NReg) = 44 + 7 mG and B,|(SReg) = - 29 + 21 mG. 



act value depends on the relative contribution of each hyper- 
fine component of the H2O 6^6—52,3 rotational maser transition. 
We adopt t he value Af^p = 0.018, which is the typical value 
reported by Vlemmings et alJ (l2002h . B\\[Gcmss] is the projected 
magnetic field strength along the line of sight, and Av^ is the full- 
width half- maximum of the I spec trum. Although the non-LTE 
analysis in IvTemmings et al. (2002) has shown that the circu- 
lar polarization spectra are not necessarily strictly proportional 
to dl/dv, using Af-p determined by a non-LTE fit introduces a 
fractional error of less than ~20% when using Eq. 1 . To increase 
the signal-to-noise ratio of the / and V spectra, we smoothed the 
data by applying a running average over three consecutive chan- 
nels. We note that the same result, albeit at lower signal-to-noise 
ratio, was obtained using the non-smoothed spectra. 

In Figure|4]we show the spectra from Stokes / and V, with al 
subtracted. Over-plotted on the V spectrum (red), we show the 
model curve that best fits the spectrum (blue). From this model 
fit, we derive a magnetic field strength of By(NReg) = 44 ± 7 mG 
and By(SReg) = - 29 ± 21 mG. The reported errors are based 
on the single channel rms using Eq.l, but similar errors were 
found using a chi-square analysis of the model circular polariza- 
tion fitted to the smoothed spectra. We also further confirmed the 
result for both maser features by performing the chi-square anal- 
ysis on spectra smoothed to different spectral resolutions (up to 
~ 0.2 km s~'). This leads us to conclude that the detection in the 
SReg is marginally significant. The flux densities of the other 28 
features, however, where not sufficient for a detection of circular 
polarization. 



4. Discussion 



_2 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 4Q 5 

1 0-1-2 
Right. Ascension offset (mas) 

Fig. 3. Maser spots of the features S.05, S.06 and S.07. The tri- 
angle sizes are scaled to the intensity of each spot emission; the 
color scale shows how the velocity varies at each spot; and the 
black lines show the direction of linear polarization for those 
spots that have survived to the polarization threshold cut. The 
size of the black lines are scaled to the linearly-polarized inten- 
sity measured for each spot. 



4.1. Kinematics 

Assuming that the maser features observed in this work, at the 
epoch of March 1" 2009, correspond to approximately the same 
features detected by D07 on November 24''' 2002, we can es- 
timate the apparent proper motion between both observational 
epochs. We compared the distance between the mean position 
of all features located in region 2. A to the mean position of all 
features located in region 2.D in our data (Figure|2]l, against the 
distance between the two corresponding points from D07. We 
present these mean positions and the offsets between the two ob- 
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Table 2. Coordinates of the center (mean) points of 
the features located in regions 2. A and 2.D, from 
Figure |2j and their offsets between the observations 
from D07 and ours. 



Reference 


Region 


Ace (mas) 


AS (mas) 


Present work 


NReg 


+2.3 


+29.2 


Present work 


SReg 


+9.7 


-44.0 


D07 


NReg 


+2.6 


+26.2 


D07 


SReg 


+2.4 


-34.6 


Offset 


NReg 


-0.3 


+2.9 


Offset 


SReg 


+7.4 


-9.4 


Offset 


Total 


7.7 


12.3 



servational epochs in Table [2] The offsets^lin a and 5 give a total 
offset of 14.4 mas which, at a distance of 1540 pc, corresponds 
to 22.2 AU. Accordingly, for the time elapsed between the obser- 
vations (2288 days), the velocity of the separation between the 
two mean points on the sky is computed to be 2.3 mas/year, or 
16 km/s. Assuming t he inclination of t he bipolar flow is /=36° 
dKastner et alJI 1992k IShure et al.l [l995). the real separation ve- 
locity is 21 + 11 km/s. 

Furthermore, we considered the overall distribution of the 
maser spots for all features that we observed. In Figure 1 .a we 
show that the spots in the NReg appear to align well with the 
direction of the nebular jet. The same property is found for the 
spots from the SReg, as can be seen more clearly in Figure [3] 
In particular, the spots from feature S.07 appear to be closely 
aligned with the direction of the jet. However, despite being co- 
linear with the jet, the H2O maser outflow velocity is signif- 
icantly lower than that r eported for the jet, vj et ~ 330 km s _1 
dSanchez Contreraslll997l) . Hence, our interpretation is that the 
masers are being entrained by the jet but are not located in the 
jet itself. As a consequence, the maser spots have a velocity gra- 
dient; with higher velocity further from the origin. The nature of 
the H2O masers in OH23 1 .8+4.2 is thus very different from that 
of the water fountain sources wher e the H?Q maser s lie at the 
tip of their fast bipolar outflows (e.g llmai et alj |2002). 

It is important to highlight, however, that a more precise 
kinematics analysis will require more accurate astrometric ob- 
servations. 

4.2. Polarization 

We detected linear polarization for three H2O maser features, all 
of them located in the SReg (S.05, S.06 and S.07). According 
to maser polarization theory, the direction of the linear polariza- 
tion can be either parallel or perpendicular to the direction of the 
magnetic field lines. It is parallel when the angle between the 
magnetic field and the direction of maser propagation (8) is less 
than the van Vleck angle ( ~55°), and it is perpendicular when 
6 is greater than this angle (Goldr eichll 19731) . The percentage of 
linear polarization is affected by and by the degree of satura- 
tion. Based on the level of linear polarization detected (~ 1%) 
we cannot conclude definitively in which regime - parallel or 
perpendicular - the masers of OH231.8+4.2 occur. 

The linear polarization angles of the features S.05, S.06 and 
S.07 display a significant scatter (Fig[T]d). The cause of this scat- 
ter is not known, but we propose three different scenarios: (i) the 
scatter could be caused by turbulence, or (ii) in the case of a 

2 These offsets occur in tangent plane project in the sky. So, strictly, 
the a and 6 listed in Table[2]correspond to the coordinate angles of these 
projections. 



toroidal magnetic field, the masers could be located at its tan- 
gent points; both would explain the different local directions of 
the field, or (iii) the scat ter could be caused by internal Faraday 
rotation (Faradav 1845). Under the effect of Faraday rotation, 
the linear polarization angle of the emission (O) is rotated by 

<D[°] = 2.02 x 10~ 2 L[AU] n e [c m - } ] B^mG] v~ 2 [Gffz] (2) 

where L is the path-length of the maser emission through a 
medium with magnetic field By and electron density n e at a fre- 
quency v. By assuming typical values L~10 AU and Z?n~50 mG- 
for the H2O masers, Eq. 2 becomes: (D~0.02 o x n e . For the 
Faraday rotation angle to be of the order of the observed scat- 
ter, an electron density of order 10 3 is required. For example, 
if n e ~2000 cm 3 , then O~40°. For an electron density of this 
order at this distance to the central star(s), the fractional ioniza- 
tion would be significantly higher then for a regular evolved star. 
However, considering the pre-PNe nature of OH231.8+4.2 this 
cannot be completely ruled out. Still, we consider either scenar- 
ios (i) and (ii) more likely than the internal Faraday rotation ef- 
fect. We note, however , that both scenari os (i) and (ii) differ from 
the results reported by lEtoka et alj d2009l) for the OH maser re- 
gion. The OH masers indicate a uniform field that seems to be 
flaring out in the same direction of the jet. However, the H2O 
and OH masers occur in very different regions. While the cur- 
rent work finds the H2O masers to be located at ~40 AU from 
the CS, the OH masers are located on a torus w ith a radius of 
~2 arcsec dZiilstra et al.ll200U lEtoka et alj|2008h . At 1540 pc 
this correspond to a distance of ~3100 AU. 

We detected circular polarization for two maser features. 
Based on these detections, we determined the strength of the 
magnetic field in each feature to be (By(NReg) = 44 + 7 mG, 
B||(SReg) = - 29 + 21 mG). If, as argued in Sect. 3.3 that the 
detection in the SReg is marginally significant, it is evident that 
the magnetic field in the SReg has an opposite sign compared to 
the field in the NReg - the positive sign indicates that the direc- 
tion of the magnetic field is away from the observer, while the 
negative sign corresponds to a direction towards the observer. 

A precise determination of the morphology of the magnetic 
field is not possible from the current work, as different field mor- 
phologies can fit the circular and linear polarization results ob- 
tained here. To reach a final conclusion as to the morphology 
of the magnetic field of OH231.8+4.2 more sensitive polariza- 
tion information is necessary. If, however, we assume a toroidal 
magnetic field (B oc it is possible to extrapolate the field 
strength to the stellar surface. By taking the separation between 
the NReg and SReg, and assuming that the CS is located cen- 
trally between the two, each H2O maser region is at ~40 AU 
from the star. Therefore, from the measured By for these regions, 
the magnetic field strength on the stellar surface (taken to have a 
radius of ~ 1 AU) must be ~ 2.5 G.This is consistent with the re- 
sult found for W43A (1.5 G: IVlemmings et al.ll2006l) . This value 
presents a lower limit if the magnetic field vs. radius relation has 
a steeper than toroidal dependence. 

5. Conclusions 

We detected 30 H 2 maser features towards OH23 1.8+4.2. 
By comparing the cu rrent work with the prior observations of 
Des murs et al.l d2007l) . we can conclude that the features are 
moving away from the central star. The average separation ve- 
locity of the masers is 21 + 11 km/s. Furthermore, the masers 

3 Note that solar-type and dipole fields have r" 1 and r~ 3 dependen- 
cies, respectively. 
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appear to be dragged in the direction of the collimated outflow. 
This likely indicates that the masers arise in the turbulent mate- 
rial that is entrained by the jet. 

We detected linear polarization for three H2O maser features. 
The large scatter in the directions of the linear polarization could 
be caused by turbulence, or could be due to the masers being lo- 
cated close to the tangent points of a toroidal magnetic field. The 
possibility of Faraday rotation has also been investigated but, un- 
less the electron density is exceptionally high (> 2000 cm -3 ), is 
ruled out. 

Based on the detection of circular polarization for two maser 
features, we determined that the strength of the source magnetic 
field is Bn(NReg) = 44 + 7 mG and B,|(SReg) = - 29 + 21 mG. 
The exact morphology of the field in the H2O maser region could 
however not be determined. Assuming a toroidal magnetic field 
(B oc r _1 ), the extrapolated magnetic field strength on the stellar 
surface is ~ 2.5 G. 

Our polarizatio n detections, together with the results from 
lEtokaet alJ (120091) . make pPN OH23 1.8+4.2 the first evolved 
star that is known to be a binary and in which the presence of 
a magnetic field is confirmed throughout the circumstellar enve- 
lope. A more comprehensive model of the magnetic field mor- 
phology and its potential evolution based on a comparison of the 
inner and outer wind will require more sensitive observations. 
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